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Abstract. In standard perturbation approaches and A-body simulations inhomo- 
geneitics arc described to evolve on a predefined background cosmology, commonly 
taken as the homogeneous-isotropic solutions of Einstein's field equations (FLRW cos- 
mologies). In order to make physical sense, this background cosmology must provide a 
reasonable description of the effective, i.e. spatially averaged, evolution of structure in- 
homogencities also in the nonlinear regime. Guided by the insights that (i) the average 
over an inhomogeneous distribution of matter and geometry is in general not given by 
a homogeneous solution of general relativity, and that (ii) the class of FLRW cosmolo- 
gies is not only locally but also globally gravitationally unstable in relevant cases, we 
here develop a perturbation approach that describes the evolution of inhomogencities 
on a general background being defined by the spatially averaged evolution equations. 
This physical background interacts with the formation of structures. We derive and 
discuss the resulting perturbation scheme for the matter model 'irrotational dust' in 
the Lagrangian picture, restricting attention to scalar perturbations. 



PACS numbers: 98.80.-k, 95.36,+x, 98.80. Jk, 04.20.-q, 04.20. Cv, 04.25.Nx 
1. A new approach to perturbation theory 

Perturbation theory is a key tool in cosmology to describe the formation of structures 
in the weakly nonlinear regime and to initialize the iV-body simulations of cosmic 
structures. The standard motivation to describe perturbations on a homogeneous- 
isotropic Friedmann-Lemaitre-Robertson- Walker (FLRW) background comes from the 
known gravitational instability of the latter. However, this background is supposed, 
implicitly, to also describe the average evolution of an inhomogeneous universe model in 
the nonlinear regime, and this is implemented as a construction principle in most known 
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relativistic, quasi-Newtonian and Newtonian perturbation schemes, and Newtonian N- 
body simulations [8]. We remark that such a construction is correct in Newtonian 
cosmology: if structures evolve on a Euclidean geometry, and if they are subjected to 
periodic boundary conditions on some large scale, then the average of inhomogeneities 
is given by the assumed background (see [7] for the proofs). 

Obviously, a Newtonian cosmology has to be considered as highly restrictive when 
one moves to the framework of general relativity, and the previous construction cannot 
be expected to work there because of: (i) the relevance of the spatial intrinsic curvature 
(the second derivatives of the metric may be significant even if the metric perturbations 
are small [T61I36T - I38U471H8] ) , together with (ii) the fact that inhomogeneities are coupled 
to the spatial curvature evolution [10], and finally (iii) the absence of a conservation 
law for the averaged intrinsic curvature pji|. Also, it has been recently shown [50], in 
the special class of scaling laws for the spatially averaged inhomogeneities (so-called 
backreaction terms), that FLRW backgrounds are not only locally but also globally 
unstable as a result of structure formation and accelerated expansion, when subjected 
to perturbations whose global contribution does not vanish. As soon as a homogeneous 
and isotropic self-graviting system is perturbed, the related inhomogeneities invoke a 
departure of its average from a FLRW background. This property is of great significance 
for a theory of perturbations: in the standard approach the background evolution is 
known to impact on the evolution of perturbations, but the converse effect, namely 
that the structure inhomogeneities also affect the evolution of the background, has been 
neglected thus far by construction^]. Note that such a property is expected from first 
principles: it expresses the fact that, in Einstein's theory, the formation of structures 
and the evolution of geometry are mutually (and generically) coupled. An evolution of 
structures on a predefined background is, in light of these remarks, an a priori restricted 
approach. 

Looking at perturbations in the universe, we can only apprehend their strength and 
their evolution correctly if we know with respect to which background they have to be 
considered. Previous work that has addressed the issue in Newtonian cosmology [581159] 
faces the drawback that, on some large scale, all averages are strictly free of backreaction 
from inhomogeneities due to the restriction to a non-dynamical geometry and the 
necessity of a torus architecture, as explained above. The same drawback remains in 
quasi-Newtonian relativistic perturbation schemes for the evolution of gauge- invariant 
variables on a predefined background, even if they take into account backreaction 
effects [39TH31 [5T| [52l [66] . A related discussion by Rasiinen can be found in [46] . 
For example, the assumption of periodic boundary conditions on initially flat space 
sections or the restriction of the scalar curvature to a constant curvature suppress any 
interaction effect between structures and the background. Beyond the usual FLRW 
perturbation scheme, Clarkson et al. have furnished in [22] a complete system of master 

| Or it has been claimed to be a negligible effect. However, arguments for a small backreaction rely 
on a weakly perturbed FLRW cosmology, and only address the issue within a limited framework (see, 
e.g., the review [24] for a clear presentation of this issue, and references therein). 
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equations that represent the general linear perturbations to Lemaitre-Tolman-Bondi 
(LTB) cosmologies. In [31] . Giesel et al. have given the full gauge-invariant relativistic 
perturbation theory holding for any background metric. 

The framework presented in this paper provides the needed tools to implement a 
background as the average over fluctuating fields, which is not based on the introduction 
of a predefined background and deviations thereof; it furnishes the evolution of the scalar 
parts of the deviation fields off this general background in a non-perturbative way. We 
expect from this improvement that we shall be able to explicitly imprint structure 
inhomogeneities into the background that eventually describes the expansion history of 
the universe without the need to invoke, e.g., a dark energy fundamental component. A 
further aim is to understand initial conditions for relativistic numerical simulations of 
inhomogeneities that are not restricted to vanish on average on a FLRW background. 

We proceed as follows. In section [2} we briefly recall the local and averaged evolution 
and constraint equations for the description of inhomogeneous dust cosmologies, and we 
provide a system of equations for the deviation fields off a general background. This is 
followed, in section |3l by a thorough discussion of the properties of this new deviation 
scheme. We conclude the paper in section |4] with some prospects. 

2. Perturbation scheme on a general background 

2.1. Local and averaged equations for inhomogeneous cosmologies 

Let us consider a globally hyperbolic four- dimensional manifold, endowed with some 
metric tensor 4 g. An irrotational fluid congruence, defined by a unit timelike vector 
field u, will be used to foliate the spacetime into a family of flow-orthogonal space-like 
hypersurfaces. We shall restrict ourselves, in what follows, to the case of a pressureless 
irrotational fluid (irrotational dust), T^ u = gu^u u , with g its energy density, as described 
in the Lagrangian picture^. In the canonical bases (d t , <9j) and (dt, dX % ), the four-velocity 
of the fluid assumes the form 

^ = (1,0,0,0), u M = (-1,0, 0,0), (la) 

and the line element is written 

ds 2 = tg^dX^dX" = -dt 2 + h ij dX i dX j , (lb) 

X 1 being the Lagrangian spatial coordinates (coordinates comoving with the fluid), and 
h the inhomogeneous three-metric of the t-constant hypersurfaces. The foliation of 
Einstein's equations with respect to u implies the well-known Raychaudhuri equation 
and Hamilton constraint: 

+ ^0 2 = -AnGg - 2a 2 + A , (2a) 

§ Greek indices refer to spacetime components, they run in {0, 1, 2, 3}, and Latin indices denote space 
components, running in {1, 2, 3}. For a presentation of the (3+1) -splitting of Einstein's equations, see 

e.g. pQGllGIl]. 
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Throughout our study, the overdot will stand for the covariant derivative (here identical 
to the partial time derivative d t ). 6 := h^Qij := ^h^hij and a := (lo^'cy) 1 ' 2 are 
the local expansion and shear rates, respectively, while TZ is the local three-Ricci scalar 
curvature of the hypersurfaces, and A is the cosmological constant that we carry along 
for the sake of generality. Those relations may be supplemented by the fluid continuity 
equation and a balance relation between TZ and a: 

g + eg = 0, (2c) 

n+-QTZ = 2{a 2 y+AQa 2 , (2d) 
3 

where the latter is obtained upon requiring equation (|2bl to be an integral of (I2al) 
(we shall simply call, hereafter, 'integrability condition' the result of this procedure). 
The local system (T5]) will be used, later on, to derive the evolution of the deviation 
fields inside any spatial domain. For notational ease, we have omitted the time and 
space dependences, but the reader should bear in mind that all variables at stake are 
inhomogeneous . 

Note : For a locally isotropic cosmology, a = 0, the system (T5]) becomes homogeneous 
(Schur-Trumper's theorem [551160] ). We recover a FLRW cosmology. In particular, 
writing the expansion rate as Oh ='■ Sdn/an with an the scale factor, the integrability 
condition ( I2dl) implies that the spatial Ricci scalar curvature follows the evolution law: 
TZh = Qk/a 2 H , where k is a constant of integration. 

Let us now consider a scalar field ip. Its spatial average performed on some compact 
domain P, contained within the hypersurfaces and transported along the fluid flow lines 
(Lagrangian averaging), is defined as (see [9] for details): 

Wv-=TT I VV det h ij( i 3 X, (3a) 

VV Jv 

with Vd := j v a/ det h~ij d 3 X the volume of the domain under consideration satisfying 
Vd/Vd = (@)x>- We shall also make frequent use of the commutation rule between 
spatial averaging and differentiation with respect to time: 

mv - &)v = (&i>) v - m v (i>) v , (3b) 

where the right-hand side reduces to zero for a homogeneous domain. Equipped with 
these relations, we can provide the Lagrangian averaging on D of Raychaudhuri's 
equation and Hamilton's constraint!]]]: 

+ ~ = -4vrG (q) v + Qv + A, (4a) 
1 -(e)l = 8nG( g ) v - { ^-^ + A, (4b) 

|| For comprehensive reviews on averaged inhomogeneous cosmologies in general relativity, we 
recommend the reading of, e.g., [H[T3[TS1|231I21|271I151I1B1|B2] (and references therein). 
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and we obtain the conservation law for the total rest mass within D and the integrability 
condition for the averaged variables as 

(QY-D + (Q)v{e)v = 0, ( 4c ) 

Wv + \ {K) v + Qv + 2 (Q) v Q v = , (4d) 
Qx> being the kinematical backreaction, 

Qv := I ((e-(e) v )%- 2 {a 2 ) v . (5) 

Equations (jlj), the averaged counterpart of d2J), will also be used in the following to 
obtain the evolution of the deviation fields in the interior of D. Qx> determines how 
the fluid inhomogeneities inside the domain globally contribute to the evolution of its 
background (equations (]4aj) and f!4bj) ). and this variable is dynamically coupled to the 
averaged scalar curvature (equation (jldl) ). Equation (l4d|) also shows that the averaged 
curvature does not individually obey a conservation law like the fluid density; rather 
a combined expression of intrinsic and extrinsic curvature invariants is conserve cfl It 
is important to notice that the background is scale-dependent: for another domain 
we have in general a different background, due to the unconstraint distribution of 
inhomogeneities. There so exists a deep correlation between the background of any 
domain and the inhomogeneities inside. This feature is habitually absent in the usual 
cosmological perturbation schemes on a global scale where perturbations are assumed 
to average out on a predefined background. 

Note : If we continuously shrink the compact domain to a point (null-homotopy) , 
((9 - (6)^)% -> and Q v -> -2a 2 . The system gj) then reduces to ©. 

Finally, it is also convenient for later discussion to introduce two of the scalar 
invariants of the expansion tensor l := cr ^ + ^@h l j, its trace and the dispersion of its 
non-diagonal components: 

I := tr(ey = 6 , II :=| (tr 2 (6^) - tr(e^)) = ^6 2 - a 2 . (6) 
The systems (J2J) and (j4j) assume the same form with these variables [9]: 

i + I 2 = 211 — AuGq + A , II = 8irGg - y + A , (7a) 

U + hn+ (211 - ^I 2 ) ' + 21 (211 - ^I 2 ) = , (7b) 
for the local one, and 

(I) v + (l) 2 v = 2 (ll) v - AttG ( e ) v + A , (ll) v = 8nG (g) v - ^ + A , (7c) 
Wv + \ (Ob + (2 (H) v ~ \ (1)1 ) ' + 2 (l) v (2 (II>„ - I (1)1 ) = , (7d) 

% Qv can be written in terms of the extrinsic curvature of the hypcrsurfaces, Kij =: —Qij, as 
Qv = (K 2 - K^K^T, - § (K%. 
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for the averaged one. We have reformulated the kinematical backreaction as Qx> = 
2 (11)^, — I (I)p to obtain the second set of expressions. This statement obviously holds 
true for the local and averaged continuity equations. 

2.2. Equations for the deviation fields 

We now provide the whole set of equations for the deviation fields off the background of 
a comoving dust domain. As is customary, we shall designate the deviation (or peculiar) 
field of any scalar field ip from its background value by 5ip := ip — (ip)x>- One would 
prefer to write 5x>ip in order to make explicit the scale-dependence of the deviations; 
however, we drop this index for notational ease. Note finally that, being a scalar 
(refer to [30] for a proof), the deviation Sip is also a scalar. 

In this paragraph, we only add a few remarks about each proposition. A thorough 
discussion follows in section [3j 

2.2.1. Deviations in density Using the local and averaged conservation laws (l2cl) and 
(J3cJ), we find the following continuity equation for the fluid density deviations, which we 
formulate in the form of a first proposition. 

Proposition la: The evolution equation for the density deviations on a compact 
domain T> is given by: 

(5g)- + (Q) v Sg = -5Q ((g) v + Sg) , (8a) 

or, equivalently, in terms of scalar invariants: 

(5g)- + (l) v 5g = -5l ((g) v + 5g) . (8b) 

Remarks : The density deviation field does not obey a conservation law like the local 
and averaged densities. We are faced with a source term involving the deviation of the 
expansion rate from its background value. By making use of the commutation rule, the 
average on V of these equations results in identities. 

For later discussion we shall prefer to use an alternative form of this proposition. 
Consider to this end the scale-dependent contrast density Ax? := Sg/g (—00 < A^ < 1), 
which is more adapted to the Lagrangian picture and the nonlinear situation (3JH], 
than the conventional definition used in Eulerian perturbation theory, 5x> '■= Sg/ {g) v 
(restricted to the same domain). By means of the local and averaged Raychaudhuri 
equations ( |2a|) and ( )4a|) . we end up with the following evolution equations for Ax>- 

Proposition lb: The evolution equations for the contrast density on a compact domain 
T> are written: 

A P = 50 (A v - 1) , (8c) 
A© + I (e) v A v - 4ttG (g) v A v = 5Q (A v - 1) , (8d) 
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or, equivalently, in terms of scalar invariants: 

A v = 51 (A v - 1) , (8e) 



A v + | (l) v A v - 4ttG (g) v A v = 2 (sil - | (l) v 51 ) (A v - 1) . (8f ) 

Remarks : We here have introduced the local contribution of fluid inhomogeneities within 
the domain, Q := |(G - (<d) v ) 2 -2a 2 = 211- § (l) v (21- (l) v ). By construction we have 
(Q) v = Qv-i and 5Q stands for the deviation of the kinematical backreaction. Taking 
the averages of these relations and using the definition of A-p and the commutation rule, 
we get identities, as it should be for a proper definition of deviation fields. 

2.2.2. Deviations in kinematical variables Using the local and averaged systems, 
( (I2ap .(l2b ]) ) and ( fj4"aj) . (14bj) ) respectively, we find the following equations for the 
kinematical deviations. 

Proposition 2: The evolution and constraint equations for the kinematical deviations 
on a compact domain T> read: 

(5Q)- + (5G) 2 + - (Q) v 5Q = -4nG5g + 5Q , (9a) 
3 

\ (Q) v 5Q = 8nG5g - Un - l -5Q , (9b) 

or, equivalently, in terms of scalar invariants: 

(51)- + (5l) 2 + 2 (1} V 51 = 2 511 - 4irG5g , (9c) 

511 = 8nG5g - -5TZ . (9d) 

Remarks : Contrary to the relation (I9dj) . the shape of fl9c|) is not identical to that of its 
local and averaged counterparts, ( 17a|) and ( 17cj) . The nonlinear character of the latter 
makes appear the extra term 2 (I)^ 51. Taking the averages of (pel) and (j9d|) . and using 
the commutation rule for the first one, we end up with identities. 



2.2.3. Integrability condition Finally, demanding equation (I9bl) to be an integral of 
( )9a|) . we get the integrability condition between the kinematical deviations and the 
intrinsic curvature, which we formulate in the form of a last proposition. 

Proposition 3: The integrability condition on a compact domain T> reads: 

(5K)- + ^(Q) v 51l+(5Q)- + 2(Q) V 5Q = -5Q {^{U) V + 2Q V + 5U + 5Q^j , (10a) 

or, equivalently, in terms of scalar invariants: 

{Sny + (l) v 51l + 2(51iy +2{l) v 511 = -1QitG51 {{g) v + 5g) . (10b) 

Remarks : The curvature and backreaction deviations are coupled through this 
integrability condition. Here again a source term involving the deviation of the 
expansion rate is present; in terms of scalar invariants, this term is the same as that of 
proposition la. The average of these relations also results in identities. 



Relativistic cosmological perturbation scheme on a general background 



8 



3. Properties of the deviation scheme 

3.1. Discussion 

We have generalized the usual dust perturbation scheme off a predefined background 
to deviations off a general background, that was obtained through the spatial average 
of inhomogeneous fields on a generic domain. The background of a pressureless self- 
gravitating system is then not restricted to follow a predefined evolution, but rather 
an evolution depending on the inhomogeneous distribution of matter and geometry. 
As we shall see further below, deviations off FLRW backgrounds are recovered for 
globally isotropic domains having a vanishing kinematical backreaction; they constitute 
a subclass of solutions within the present framework. Here, we prefer to speak in 
terms of deviations rather than perturbations since, by construction, this scheme is non- 
perturbative (no approximation or linearization have been performed). A consequence 
of this property, and accordingly a second interesting feature of this approach, is that 
we do not get any constraint on the strength of deviations, apart from the requirement 
of regularity of the solutions. 

It is also worth noting that this scheme only functionally depends on a metric 
(via the domain of averaging); all the equations outlined above keep the same form for 
any spatial metric. The role of the three-metric was indeed entirely implicit for the 
derivation of our scheme: we did not need to compute an averaged metric from the 
local one in order to eventually obtain the dynamics of the deviation fields off a general 
background. This is a nice feature, and one can choose any spatial metric and end up 
with a solution. We recognize that the disadvantage is to be able to deal only with 
the scalar modes of the deviations, and not the vectorial and tensorial ones. However, 
we think that the results expounded here clearly constitute a first useful step for an 
understanding of cosmological deviations off a general background. 

The generalization we have proposed shows that the kinematical backreaction, 
which encodes the global contribution of fluid inhomogeneities within a spatial domain, 
impacts in several ways on the dynamics of the deviations: not only Qd affects the 
evolution of deviations through that of the background, as it is implied by the coefficients 
like (O)^, in the left-hand side of equations ( l8dl) . ( |9al) and (l9b|) . but it also acts as a 
source through the peculiar-backreaction SQ (right-hand side of fl8dl) . fl9aj) and f )9b|) ). 
Comparing the evolution of averaged inhomogeneous cosmologies to that of (scale- 
dependent) FLRW models, it appears natural that the extra term involved in the 
evolution of density deviations off a general background is precisely the kinematical 
backreaction (see subsection 13 .4H . We also stress a difference between a general domain 
and a FLRW domain: the spatial curvature of the former is in general not given by a 
constant-curvature model. The different evolution histories of the background curvature 
influence the dynamics of the background, which modifies in turn the evolution of the 
deviation fields — again in comparison to that of a FLRW background. In addition, 
inhomogeneities in geometry also act as a source in the evolution equations (the term 
571 in equation ( )9b|) V 
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3.2. Steps toward an exact resolution 

Let us first consider the systems for the local variables, the averaged variables, and the 
deviation variables individually: 

• the local system ([2]) needs one additional relation to be solved, since we deal with 
three independent expressions for four variables (g, 0, a and TV). The hierarchy 
would continue with the evolution equation for the shear, but it will never be closed 
on the level of ordinary differential equations (see [9] and also [34]): 

• the averaged system (j3J) also requires a last relation (one on the averaged variables) 
since in this situation we also have three independent expressions for four variables 
((g) v , (Q) v , Qv and (TZ) V ) (see ibid.)] 

• finally, the deviation field system (JSj) calls for three additional relations, since we 
have three independent equations for now six variables (5g, (g) v , <56, (Q) v , SQ 
and 5TV). 

However, the two last systems are obviously related to the first one. Given any closure 
relation for the local system, we can solve all sets of equations and eventually obtain an 
exact solution for the deviation scheme. One can think of choosing a specific form for 
the spatial metric, or of giving a constraint on kinematical or geometrical variables (e.g. 
considering a 'silent universe model' [2]). For instance, considering a domain endowed 
with a spatial LTB metric, we are able to construct its background and compute its 
averaged evolution. According to our scheme, deviations from this background have 
to be understood as local perturbations giving back the local LTB metric. This can 
be applied to any synchronous metric. Although maybe not suited to furnish relevant 
physical models, this first procedure may provide interesting toy models. 

If we instead specify a constraint on the averaged variables on a chosen scale, 
we cannot end up with a solution for the deviation scheme. Nevertheless, by doing 
so we reduce the needed number of extra relations to one in order to exactly resolve 
it, since there will remain four variables (Sg, 50, 5Q and 5TZ) for the three-equation 
system (jSJ). The second approach to get a solution may then be realized by reducing the 
space of possible backgrounds, e.g. with scaling laws [T2|l50]. or particular effective state 
equations P2JII9] for the backreaction terms, or multiscale partitionings combined with 
closure assumptions, e.g. [6lJ|63j|6l] > an d then by considering the resulting equations for 
the deviation fields. Although being straightforward to get working models of structure 
formation, this latter approach shall always call for physical verification of the closure 
relations used. 

3.3. Definition of a global physical background and deviations thereof 

We dedicate this subsection to the reformulation of the deviation field system (jSJ) with 
the help of a spatial metric comoving with the global physical background, that we shall 
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3.3.1. The global physical background in cosmology We have discussed thus far the 
evolution and constraint equations for the deviation fields off a general physical 
background, obtained from the spatial averaging procedure. Let us now consider a 
compact spatial domain of the universe — we shall call it E — that we may assume to 
cover the homogeneity scale, namely that spatial scale beyond which all averages do no 
longer depend on scale. Such an assumption, common to cosmology, is not necessary, 
but it enables us to define a scale over which we think the universe is representative (any 
larger scale would not provide new insights) and to have a more transparent frame of 
comparison with the standard cosmological model and its usual perturbation schemes. 
We emphasize, again, that the following reformulation does not rely on the existence 
of a homogeneity scale, and it may equivalently be employed to describe globally 
inhomogeneous models, in which case E would cover the whole spatial manifold, e.g. a 
spherical space without boundary (see [11] for the average properties of such universe 
models). 

We shall identify the scale of homogeneity as the one where the global physical 
background is defined (in practice such a scale would correspond to statistical 
homogeneity — see [57] for an examination of the subject). The idea realized here 
corresponds to the ABS ("Average Background Solution") discussed by Kolb and 
collaborators [3TII3H] - All averages indexed by E shall then refer to 'global' averages 
that define our background. 

3.3.2. A globally volume-preserving metric Let us introduce, on E, the conformally 
rescaled Riemannian three-metric h as 

h i:j := 4 hij , h lj := a E 2 h ij , = = 5) , (11a) 

where is the dimensionless effective scale factor of the global domain E, defined as, 
and satisfying: 

-=(£f' 3 !^ (e> - (llb) 

with Ej the domain at initial time. The definition (II lap guarantees that the metric h 
conserves the volume of the domain E: Vz := J E (det h^) 1 / 2 d 3 X = a E 3 Vs- We may say 
that h stands for the spatial metric "comoving" with the global background. Remark 

+ The formalism propounded hereafter may as well be viewed as a reformulation of our scheme for 
spatial coordinates rescaled by the scale factor of the global background. Although reminiscent of 
the standard procedure of introducing 'comoving coordinates', we shall not pursue this possibility 
since our approach is coordinate-invariant in the hypersurfaces, and such coordinate changes would 
add no physical insight here. Comoving coordinates make sense if a global coordinate system, e.g., 
on a constant-curvature domain, can be introduced. In general this is not possible, and a conformal 
transformation of local coordinates seems unnecessary (it may be of technical help in calculations). 
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that, the conformal three- metric being inhomogeneous, we still have local variations of 
volume elements in h. 

We now propose to construct, in this metric, all the scalar fields we shall need for 
the reformulation of our scheme. Upon defining the corresponding expansion tensor as 

ey.= ^h ik h kj , (12a) 
we are able to write, with the help of expressions (II lap and (lllbj) : 



e i = 9 i - I < >e h's = o\ + 1(0 - (0} s ) h'j . (12b) 
The trace of this equality and the average of the resulting expression, respectively, yield: 

e = e-(0) s = 5e, (e> E = o. (12c) 

The first relation reveals that pinpoints the deviation of the local expansion rate from 
the h-background expansion rate of the global domain, and it thus defines the peculiar- 
expansion rate of the latter. These expressions are consistent with the stationarity of 
the h-background of E, (0}s = 0, and with the existence, in general, of local variations 
of volume elements evaluated with h, 9 ^ (6 vanishes only if the global domain is 
homogeneous) . 

By means of equation (112b|) . the traceless part of the 'tilde' expansion tensor is the 



'tilde' shear tensor of the fluid defined as 5* ;= cr l which implies a 2 = a 2 . Concerning 
the spatial curvature of the global domain, a straightforward calculation, calling for the 
use of equation f ilial) , results in 71 = a^lZ. 

Finally, the tilde energy density g is obtained by considering the fluid conservation law 

^+0^=0, (13a) 

which is simply the counterpart of ( l2c|) in the spatial geometry generated from the 
conformally rescaled metric. Writing 

Q , (Q)h ( Q Y (£)s 



= _(0) _£ + A£^ = _ * (13b) 

where we have used the local and averaged conservation laws for the second equality, 
we get 

' gyf - = 0. (13C) 



Multiplying both sides by the initial averaged density, we conclude 

g := ^J- (g)^. , (f?)s = (g) s . = const. (13d) 

The last relation is naturally expected: the total rest mass within the global domain 
being conserved, the matter tilde density has to remain globally constant in the frozen 
volume Vs- (Consider the invariant total rest mass M. — M. — Vz (g)^; its conservation 
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indeed implies = 0.) We also remark that o is the only scalar field we shall use that 
is not affected by the conformal rescaling. This characteristics is due to the fact that h 
is built such to provide a volume-preserving metric and thus only transforms quantities 
related to the trace part of tensors. 

In terms of scalar invariants, we are able to write the following equalities: 

I = (I) s +T, II = 5 + \ (I) 2 S + | (I) E T, (14a) 

where we have introduced the scalar invariants of ©^ as I := tr(0*.) = and 
II := ~(tr 2 (0 l ) — t^© 1 -©^)) = |0 2 — cr 2 . Averaging over £ the last relation, and 
inserting the result back into it, gives^ 

II = (II> S + H - <H) S + | <I) S T. (14b) 

I = 51 represents the deviation field of the first local scalar invariant from its 
average on the global domain, and II — (II) s + § (Ogl = is that of the second 
local scalar invariant. Using these relations, we also recast the local contribution 
of fluid inhomogeneities into Q = |© 2 — 2a 2 = 211 = Q, and the global one into 
Qe = |(© 2 )e - 2(a 2 )s = 2(fl) s = <2e (therefore 5Q = 5Q = 511). 

3.3.3. Deviations off' a global physical background From the propositions given in 
section [2] a set of corollaries follows. It determines the evolution and constraint equations 
for the deviation fields off the global background, as expressed in a globally volume- 
preserving metric. From this point of view, deviations do not 'see' and are not affected 
by the volume deformation of the background. We write the system in terms of the tilde 
scalar invariants, which can be straightforwardly found by using the above definitions. 

Corollary 1: The evolution equations for the contrast density on a global background 
£ read: 

A S =T(A S -1), (15a) 

A E + 2 — A 2 - 4vrG^A E = 2 511 (A E - 1) . (15b) 
as a s 



Remark : We have used the equality As = 5g/g = 5q/q = As and the definition (lllbj) . 



Corollary 2: The evolution and constraint equations for the tilde scalar invariants on 
a global background S are written: 

I+T 2 + 2— T=2 5fl-47rG^f , (16a) 

(Zs C^s 

55 + 2^1 = STrGf -- 5 4- (16b) 
as a s 2 a s 

* To be rigorous we should write lipl^ the average over E of any tilde scalar field, with Ms := 
(l/Vs) J s . \i det hij d 3 X . However, one can easily check that, for any spacetime scalar field, i.fe = (.) s . 



Relativistic cosmological perturbation scheme on a general background 13 
Corollary 3: The integrability condition on a global background E is given by: 

+ 2(jH) . + 12 * ai = -T f |% 4- 4<fl) E + ™ + 2 «A . (17) 
a E as \3 a s a E y 

5.^. T7ie limit of a FLRW background 

FLRW cosmologies are recovered only if, for any compact region T> lying in the interior 
of E, we ask for the vanishing of Qv (the integrability condition (14d|) then imposes 
1Z to be the Friedmannian curvature). In this situation, any scalar field equals its 
background value; the local and averaged systems (J2J) and (j3J) are identical, and there 
do not exist deviations over E. The existence of deviations, and hence the possibility 
to form structures, then demands the global domain not to remain locally isotropic. 
In other words, we need to abandon the strong cosmological principle in favor of, for 
instance, a weaker version that defines a scale of homogeneity (see subsection I3.3p . We 
may for example consider E to follow globally, and not locally, a FLRW evolution, 
requiring the kinematical backreaction to vanish on the homogeneity scale E, and 
therefore on any larger scale. Such an assumption ensures that the background is 
globally Friedmannian, but at the same time does not prevent local inhomogeneities 
to live inside the global domain. In this picture, the cancellation of <2s would be the 
result of an exact compensation between the expansion variance and the average of the 
shear squarec^]. 

Facing such a global domain constrains the averaged system with the condition 
Qy, = 2(II)£ = 0, and allows us to write the second-order differential equation of the 
contrast density ( 115bl) as follows: 

A s + 2^A E - 4ttG%A s = 2I~I (A s - 1) , (18a) 
as a s 

where as is the scale factor of the background satisfying: 

3^ = -4ttG% + A , 3 V = 8ttG% - ± <% + A . (18b) 
as a s V a s/ <4 2 a s 

We recover with (118j) the standard framework of density deviations off a FLRW 
background cosmology j6], which therefore constitutes a subclass of solutions of our 
deviation scheme. We leave it to the reader to simplify the other deviation field equations 
(corollaries 2 and 3) using the above condition (511 = II). Let us further suppose that 
the fluid inhomogeneities weakly contribute to the local kinematics of the deviation 
fields: Q = Q ~ 0. We can then neglect the quadratic invariant II and end up with: 

As + 2— As - 4ttG^ = , (19a) 
a s a s 

which gives the evolution of the first-order Lagrangian (relativistic) density 

perturbations off a FLRW background. The same expression is obtained for the 

jj This is what happens, for instance, for spatially averaged zero-curvature LTB models. 
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evolution of the linear density perturbations in standard perturbation theory. This 
equation corresponds to the linearization in As of (]15b|) . and it is solved by the 
exact relativistic form of Zel'dovich's approximation [33j|65], systematically derived in 
a forthcoming paper |20j . 

We finally take advantage of this small- deviation picture to make a digression about 
the usefulness of the contrast density A s over the density contrast 5 S = A s /(1 — A s ). 
Expressing relation (11 9 aft in terms of 5% yields: 

5 E + 2^ E -47rG(£) E 5 s 

+ (fe + 2^5 S - 87rG (e) s 5 S - 47rG ( e ) s 51) - 2tf| = ■ (19b) 

This illustrates that the solution to the linear equation (I19a[) for As substantially goes 
beyond that for <5s, 5 S + 2(ds/as)<5s — 47rG 5 s = 0. Hence, first-order Lagrangian 
(relativistic) deviations off a FLRW background already involve nonlinearities in the 
dependent variable 5s, which demonstrates the inherently nonlinear character of a 
Lagrangian perturbation approach. Remark at last that the density contrast coincides 
by construction with the density deviations 5g evaluated in the globally comoving metric; 
it is simply the relativistic extension of the density deviation field used in the standard 
Eulerian cosmological perturbation theory [44] (see [H1E1I26] for other remarks on the 
Lagrangian picture versus the Eulerian one in Newtonian cosmology). 

4. Concluding remarks and outlook 

We have generalized, in the present paper, the dust scalar perturbation scheme off 
a predefined background to deviations off a general background. The kinematical 
backreaction, which determines the global contribution of fluid inhomogeneities within a 
generic domain, is at the very core of this non-perturbative scheme: it not only influences 
the dynamics of a general background, as it is well known, but it also explicitly impacts 
on that of the deviation fields. Our longterm expectation from this improvement is 
to be able to describe large-scale structure formation uniquely from the existence of 
inhomogeneities and without the need to invoke a dark energy fundamental component. 

Our scheme may be exactly solved either by considering a specific three-metric, 
or by imposing local dynamical constraints, or by restricting attention to subclasses of 
backgrounds and then constraining the deviation fields. Another, and we think the most 
promising, strategy to solve the deviation equations would be to develop an iterative 
procedure. The reason why an iterative procedure takes better care, compared to a 
perturbative approach, of the nonlinear character of the proposed scheme is obvious: a 
perturbation point of view runs into contradiction due to the fact that the background 
(the zeroth-order solution) is generally modified by the kinematical backreaction (a 
second-order term). Hence, for situations where the backreaction term does not vanish, 
the notion of, e.g., first-order deviations off a general background would be ill-defined. 
An iterative point of view also entails methods that are known to numerical simulations, 
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and we expect that the simplest application of the presented scheme is numerical in 
nature. This relativistic Lagrangian procedure would consist, for the O th -level, of taking 
the Zel'dovich approximation (11 9 aft and computing the l st -level backreaction term. This 
first implementation has been depicted in [18] (see section 7.3). Reinjecting this 1 st - 
level solution, at each step of time, into the background, perturbating this latter and 
solving the equations of the above corollaries, would then drive the 2 nd -level kinematical 
backreaction, and so on for the ensuing levels of iteration. This process should clearly be 
carried along beyond the first levels of iteration, in order that the inherently nonlinear 
character of structure formation, and its effects, shall be taken into account. Another 
promising approach would be to use this procedure in the framework of the gradient 
expansion treatment of inhomogeneities. Developed in [25l[53] and used in [35] to study 
backreaction, this technique — contrary to perturbation theory — seems to take better 
care of the impact of the small-scale nonlinear effects on the larger ones. A recent 
study [29] has shown that, already for gradient expansion quantities of the 4 th - level, the 
effects of backreaction grows to 5-10% of the background. A realization of this iterative 
strategy is the subject of forthcoming work. 
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